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Abstract. We establish the solvability of second order divergence type par- 
abolic systems in Sobolev spaces. The leading coefficients are assumed to be 
only measurable in one spatial direction on each small parabolic cylinder with 
the spatial direction allowed to depend on the cylinder. In the other orthog- 
onal directions and the time variable the coefficients have locally small mean 
oscillations. We also obtain the corresponding VVjJ -solvability of second or- 
der elliptic systems in divergence form. Our results are new even for scalar 
equations and the proofs simplify the methods used previously in jl^ . 



1. Introduction 

This paper concerns the unique solvabihty of divergence type parabolic and 
ehiptic systems in Sobolev spaces when the leading coefHcients are in the class of 
variably partially BMO (bounded mean oscillation) functions. The parabolic system 
we consider has the form 

Vu- \u = dW g + f, (1.1) 
where A > is a constant, g = {qi, §2, - ■ ■ , gd)^ ^^id 

Vu^ -ut + Dc{A°'^Df3u) + Da{B"'u) + B°'DcU+Cu. 

The coefficients A""^ , B" , C are m x m matrices, which are bounded and 
measurable, and the leading coefficients A"^ are uniformly elliptic. Note that 

are (column) vector-valued functions defined on 

{S,T) X M'* = {{t,x) : t e (S',r), X = [xi,-- - ,Xd) G R'^} , 

where — oo < S < T < oo. For given g and / G Lp, I < p < oo, we seek a unique 
solution u in the parabolic Sobolev space Hp (for a definition of Hp, see Section ||). 
We also consider the following elliptic system 

Cu - Xu= div g + f, (1.2) 

where 

Cu = Da[A°'^Dpu) + Do,{B"u) + B"DaU+ Cu. 

In this case A"-^ , B", B°', C, g, and / are independent of t and satisfy the same 
conditions as in the parabolic case. Naturally, the solution space is Wp. 
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There are many papers concerning elliptic and parabolic equations/systems in 
Sobolev spaces with VMO (vanishing mean oscillation) or BMO type coefficients. 
Chiarenza, Frasca, and Longo first proved the interior estimate for non-divergence 
form elliptic equations with VMO coefficients |^ . Then the solvability of elliptic 
and parabolic equations in Sobolev space were presented in and j^. These are 
the earliest papers about non-divergence type equations with VMO coefficients. 
For divergence type equations with VMO coefficients, results of similar type were 
obtained in 0. Later, Byun and Wang studied in their papers (see, e.g.,[^ |j 
and references therein) divergence type equations/systems with BMO coefficients 
in non-smooth domains. Krylov also treated BMO coefficients in |2^, where he 
gave a unified approach to investigating the ip-solvability of both divergence and 
non-divergence form parabolic and elliptic equations with coefficients BMO in the 
spatial variables (and only measurable in t in the parabolic case). For other related 



results, we also refer the reader to |2a, 26, 21, n5l H, p|, and references therein 



To explain the class of coefficients in this paper, we first mention partially BMO 
coefficients, which are characterized as having no regularity assumptions with re- 
spect to one (fixed) variable and having locally small mean oscillations with respect 
to the other variables. This class of coefficients was first introduced in where 
the Wp-solvability of elliptic equations in non-divergence form were obtained by 
adapting some ideas in p^p[ Since then, non-divergence type equations /systems 
with partially VMO/BMO coefficients have been considered in |l6[ |l^, 
Partially BMO coefficients are quite general so that they include VMO coefficients 
as in |6[ 0, P as well as BMO coefficients as in ^. As to divergence type equa- 
tions, the authors of this paper proved in jl^ the W^^-solvability of elliptic equations 
with partially BMO coefficients. Then parabolic equations as well as systems in 
divergence form were treated in [p^ ^ . 

In this paper, we deal with the class of variably partially BMO coefficients, which 
is a generalization of partially BMO coefficients. This class of coefficients was first 
introduced by Krylov in ||2^ for elliptic equations in non-divergence form in the 
whole space and p G (2, oo). Variably partially BMO coefficients are measurable in 
one spatial direction and have small mean oscillation in the other directions via a 
diffeomorphism on each small cylinder (or ball in the elliptic case) . DifFeomorphisms 
may be chosen differently for each cylinder, so the direction in which coefficients 
are only measurable (have no regularity assumption) may vary from one cylinder to 
another. In other words, there is no global fixed direction, with respect to which the 
coefficients are only measurable. It is easily seen that the class of partially BMO 
coefficients is a special case of variably partially BMO coefficients with the identity 
diffeomorphism. Later, non-divergence type parabolic equations with similar type 
of coefficients and any p G (1, oo) were dealt with in pl| . 

Having variably partially BMO coefficients in hand we establish in this paper the 
corresponding results of and for divergence type equations, which generalize 
the results of Moreover, in contrast to |ri|, ^ we deal with systems as 
well, so our results extend all results in to the system case. The main steps 
are L2-oscillation estimates of the derivatives of solutions, and then applying a 
generalized Fefferman- Stein theorem proved by Krylov in ||2^. However, there are 



*In fact, the authors of jl^ considered partially VMO coefficients, which are only measurable 
in one fixed variable and have vanishing mean oscillations in the other variables. However, in the 
same spirit as in the proofs there also work for equations with partially BMO coefficients. 
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additional difficulties due to the divergence structure of the equations/systems and 
the fact that coefficients are (locally) only measurable in one direction. 

To overcome the difficulty due to the divergence structure, in jl2| we used a scal- 
ing argument. Roughly speaking, we considered a rescaled function u{^~^xi,x') 
instead of u{xi, x') to get a priori estimates of Du, where is a large constant and 
{xi,x') e M''. The coefficients considered in this paper, as noted above, have no 
specific fixed direction to which we can apply the scaling argument. This prompted 
us to develop a new method, the key step of which is to estimate the mean oscilla- 
tions of X]j=i o}^ DjU and DiU, i = 2, ■ ■ ■ ,d, instead of the full gradient of u, if the 
given equation is 

Di{a^^ Dju) = div g. 

Applying our new method to the equations in , it not only removes the necessity 
of the scaling procedure, but also simplifies the proofs there. Moreover, the method 
allows us to treat systems, whereas in we were only able to deal with scalar 

equations due to the fact that a certain change of variables had to be used. 

Comparing to elliptic equations considered in ||l^, another obstacle in the par- 
abolic case is in the estimate of || tit Hl^. In contrast to the case of non-divergence 
form equations, the estimate of ||ttt||L2 does not follows directly from those of spa- 
tial derivatives of u. To circumvent this obstacle, in Lemma 3.2 we use an iteration 
argument combined with suitably chosen weights. For equations with symmetric 
coefficient matrices, a simpler proof can be found in | [Tc| |. 

Unlike |l2|, where equations are considered in the whole space, a half space and 
a bounded domain, here we only concentrate on equations in the whole space for 
the simplicity of the presentation. For a discussion about different approaches for 
equations with VMO, BMO, or partially BMO coefficients, see and references 
therein. 

The paper is organized as follows. We introduce some notation and present the 
main results in Section 0. In Section S we prove some preliminary results, which 



are necessary in the proofs of the main results. We prove in Section |J Theorem 2.2 
The other main results are easily derived from Theorem 2.2. Finally in Section |5| 



we make some remarks on elliptic systems with some less regularity assumptions 
on diffeomorphisms. 



2. Notation and main results 

2.1. Notation and function spaces. We begin this section by introducing some 
notation, which will be used throughout the paper. Let m,d > 1 be integers. A 
typical point in R'^ is denoted by x = {xi, • • • , Xd) = {xi, x'). We set 

By Du and D^u we mean the gradient and the Hessian matrix of u. On many 
occasions we need to take these objects relative to only part of variables. In such 
cases, we use the following notation: 

Dx' — '^x' : Dxi x' — '^^yx' i Dxx' — '^xx' 7 

where, for example, Dxx'U means one of Ux^x/i, a = 1, • • • , d, /3 = 2, • • • , d, or the 
whole collection of them. 
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Throughout the paper, we always assume that 1 < p < oo unless explicitly 
specified otherwise. By N{d,m,p, ■ • • ) we mean that is a constant depending 
only on the prescribed quantities d,m^p, ■ ■ ■ . 

For a function f{t,x) in M''+^, we set {f)v to be the average of / over an open 
set V in M'^+i, i.e., 

{f)v = J^f{t,x)dxdt ^ f{t, x) dx dt, 

where is the d + 1-dimensional Lebesgue measure of T). 
For — oo < S < T < oo, we denote 

nl{{S,T) X R'^) = (1 - A)i/2w^pi^2((5,T) X R'^), 

H-i((5,r) X R'^) = (1 - Ay/^LpiiS,T) x R'^). 

For any T £ (— oo, oo], we use 

Mt = (-00, T), M^+i = Rt X R-^ 

to abbreviate, for example, ip((-oo,T) x M'') = Lp{Rp'^). When T = oo, we 
frequently use the abbreviations Lp = Lp(R'^+i), Hj, = 7i^(R'*+i), etc. 
Set 

B^x') = {ye R"-' : \x' - y'\ < r}, Br{x) = {y : \x ~ y\ < r}, 
g;(i, x) = {t- r^,t) X Bl{x'), Qrit, x)^{t-r^,t)x 

and 

s; = b;(o), Br = Br{o), Q; = g;(o,o), = q;(o,o). 

As above, |_B^|, \Br\, \Q'r\, and \Qr\ mean the volume of B'^., Br, Q'^, and Qr 
respectively. 

For a function g defined on R'^+i, we denote its (parabolic) maximal and sharp 
function, respectively, by 

Mg{t,x)= sup / \g{s,y)\dyds, 

QeQ.{t,x)eQ J Q 

g*{t,x)^ sup + \g{s,y) - {g)Q\dyds, 
QeQ:{t,x)eQ J Q 

where Q is the collection of all cylinders in R'^+i, i.e. 

Q = {Qr{t, x) : {t, x) e R'^+\r e (0, oo)} . 

2.2. Main results. We state our assumptions on the coefficients precisely. We 
assume that all the coefficients are bounded and measurable, and A"^ are uniformly 
elliptic, i.e. there exist S G (0, 1] and K > 1 such that for any vectors = (Ca) ^ 
R", a = 1, • • • , d and any {t, x) € R'^+'^ we have 

dm dm 

^EEi^"i'^ E E<(*'-)^^^^ (2-1) 

a— 1 i—1 a,P—l — l 

\A°''^it,x)\<5-\ \B°'{t,x)\<K, \B°'{t,x)\<K, \Cit,x)\<K, 

where a, (3 — 1, 2, • • • ,d. 

Denote by A the set of md x md matrix-valued measurable functions A = 
[A°'^{yi)) of one spatial variable such that ( |2.l[ ) holds with A in place of A. 
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Let vj/ be the set of C^'^ difFeomorphisms -0 : M'' — > M'' such that the mappings 

ip and (p ~ Tp^'^ satisfy 

\Di;\ + \D^i>\ < S-\ \Dc^\ + \D^(j)\<S-\ (2.2) 

Assumption 2.1 (7). There exists a positive constant Bq E (0, 1] such that, for 
any parabohc cyhnder Q of radius less than Rq, one can find an A Cz A and a 
tp = (■01, • • • , 0d) g ^ such that 



\A(t,x) - Aii;iix))\dxdt < -f\Q\. (2.3) 

'Q 

Now we state the main results of this paper. Our first theorem is about the 
solvability of in R^+^ 

Theorem 2.2. Let p e (l,oo) and T e (-oo,cx3] and u S 'A^p(]R^+^). Then there 
exist constants 7 = j{d,m,p, S) > 0, and Aq > and N > 0, depending only on d, 
m, p, Rq, S and K, such that under Assumption \2. \ (7^ the following assertions 
hold. 

(i) For any uG Ti.p^U.'lf^^), we have 

^IIHIl,(rJ+1) + ^/^P"IIlp(R^+1) + KIIh-1(R^+1) < N{^+l)\\Vu- AM||jj-i(j5.d+i) 

(2.4) 

for all A > Aq. 

(ii) For any A > Aq and f, g € Lp(M.^^), there exists a unique u € 7ip(M^"^"'^) 
solving 

Vu— Xu = div g + f 
inM.'^^. Moreover, u satisfies the estimate 

^IIHlL^(R^+i) + VA||i:'u||^^(Kd+i) + ||?/t||H-i(gd+i) < iV\/A|lg|i^^(]Rd+i)+7Vy|^^(„d+i^. 



The next result is regarding the initial value problem of ( |1 . l[) . For —oo<S< 
T < cx) we define HlpHS, T) x Vl) to be the subspace of nl^p{{S, T) x Vt) consisting 
of functions satisfying uxt>s G 7ijp((— oo,T) x f2). 

Theorem 2.3. Let p E (l,oo), T e (0,cx3). Then there exist s a constant 7 > 
depending only on d, m, p and S, such that under Assumption \2.^ (j), for any f, 
gG Lp{{0,T) X R"*), there exists a unique uE HpdOjT) x R'*) satisfying 

Vu = div g + f 

in (0, T) X R"^. Moreover, there is a constant N depending only on d, m, p, T , Rq, 
S and K such that 

ll^*llwi((0,T)xR'') < N (||/IILp((0,T)xR<') + llffllLp((0,T)xR<')) • 

Indeed, by considering v e^'^'^''+^)*tt instead of u the operator V becomes 

V — (Aq + 1)/. Now we extend / and g to be zero for t < and solve the system for 

V in R^"*"^ using Theorem 2.2. By the uniqueness, we have v — when t < 0. Thus 
u ~ e.{^o+i)t^ solves the original initial value problem and the estimate follows as 
well. 

As a consequence of Theorem we obtain the W^J -solvability of elliptic systems 
(1.2) with variably partially BMO coefficients with locally small BMO semi-norms. 
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Assumption 2.4 (7). There exists a positive constant Rq e (0,1] such that, for 
any ball B of radius less than Rq , one can find an A £ A and aip — {ipi , • • • , G ^ 
such that 

' \A{x)-A{M^))\dx<^\B\. 



B 

Theorem 2.5. Let p G (l,c»). Then there exist constants 7 — 'y{d,m,p, S) > 0, 
and Xq > 0, N > depending only on d, m, p, K, 6 and Rq such that under 
Assumption \2.^ ('y) the following assertions hold, 
(i) For any u E Wp{M.'^), we have 

A||m|Ilp(R'') +\/A||i:'u||i^(Rd) < iV(\/A + Au||jj-i(Rd) 

for all A > Aq. 

(a) For any A > Aq and f,gE Lp{S.'^), there exists a unique u £ W^iW^) solving 
( |l.2| ) in . Moreover, u satisfies the estimate 

MIuWljri) + ^\\Du\\l (jga) < iV\/A||g||L (Rd) + TVyii (Rd). 



Theorem |2.5| is deduced from Theorem \i.2\ by using the idea that solutions to 
elliptic systems can be viewed as steady state solutions to parabolic systems. We 
omit the details and refer the reader to the proof of Theorem 2.6 |2^. In Section 
^ we shall give an outline of the proof of Theorem 2.5 under a weaker regularity 
assumption on ip and cj). 

Remark 2.6. The Hp-solvability results in this paper admit the extension to the 
mixed norm spaces Ti^ p by following the idea in see also, for instance, [|8| 
and [|lO[. Here we do not pursue this, and leave it to interested readers. 



Remark 2.7. As an application of Theorem |2.2| and 2.5, one can obtain the solv- 
ability of parabolic and elliptic systems on a half space or a bounded Lipschitz 
domain with a small Lipschitz constant, with either the homogeneous Dirichlet 
boundary condition or the conormal derivative boundary condition. For systems 
on a half space, near the boundary we require A"^ to be measurable in the normal 
direction and have locally small mean oscillation in the other directions. For sys- 
tems on a Lipschitz domain, near the boundary we require A"^^ to have locally small 
mean oscillation in the spatial directions (and measurable in t in the parabolic case; 
cf. Theorem 5.1 |l3|). In both cases, A"'^ is also assumed to satisfy Assumption 



2.1 (or Assumption 2.4 in the elliptic case) in the interior of the domain. We omit 
the detail and refer interested readers to the discussions in fl^ ]. 

3. Estimates of mean oscillations 

In this section we assume B = B = and C = 0. The main objective of this 
section is to estimate the L2-oscillations of solutions to Vu = div g, which is the 
key ingredient in the proofs of our main results. We start with the well-known 



7i2-solvability of (LI) with measurable coefficients. 
Lemma 3.1. 

(i) Let T G (-00, 00] and A > 0. Assume u G TiliR'^'^^) and Vu - Xu = div g + f, 
where f,g(z L2(R^^). Then just under the uniform ellipticity condition (with no 
regularity assumption on A"^^), there exists a constant N = N{d,m,S) such that 
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If X = and f=0, we have 

(ii) For A > and any f,g(z L2{M.^^), there exists a unique u G 7^2(^7'^"'^) solving 
Vu- \u = dYvg + f in K^'^ . 

(Hi) Let T G (0,cxd). For any f,gE L2{{0,T) x M'*), there exists a unique u G 
7t!^((0,T) X R'') solving ru = div g + f m {0,T)xR'^. Moreover, 

ll'"l|-Hi((0,T)xR<') < ^I15||l2((0,T)xR'') + iV||^|2,2((o,T)xR'i), 

where N = N{d,m,S,T) > 0. Iff= 0, we have 

II^HlL2((0,T)xR'i) < ^I|5||l2((0,T)xR<')i 

where N = N{d, m,S) > 

We recall the following Caccioppoli-type inequality for parabolic systems in di- 
vergence form. 

Lemma 3.2. Let < r < R < 00. Assume u G Ti-l and Vu = div g + f in Qr, 
where f,g(z L2{Qr). Then there exists a constant N = N(d,m,6) such that 

PHIl2(Q.) < N (||5||l2(Qh) + rm\L,(Qn) + - r)-'ML,iQn)) ■ 

Proof. We provide a sketchy proof for the sake of completeness. Take a ^ G 
such that 

f 1 on Qr 



c 

and 



on R''+'-\{-R^,R^) X Br 



\DC\<N{R-ry\ \Ct\<NiR- 
After multiplying both sides of the system by C^n and integrating on Qr, we get 

/ ut ■ v/;'^ dx dt + I Da{uC'^) ■ {A°''^Df}u)dxdt = - {diY g + fj ■ uC"^ dx dt. 

JQr JQr JQr 

Integrating by parts and using Young's inequality yield 

/ D^{u)-{A'^'^Dpu)(:'^ dxdt 
JQr 

<N f {{R-r)-^\u\^ + {R-rf\f\^ + \g\^) dxdt + ^ f \Du\^C^ dx dt, (3.1) 

JQr JQr 

where N — N{d, m, 6) > 0. To finish the proof, it suffices to use and absorb 
the last term on the right-hand side of (J^^ to the left-hand side. □ 

On account of the above lemma, we make a frequent use of the following ar- 
gument. Let A"^ = A°'^{xi), that is, they are functions of G M only. Also 
let u G and assume Vu = in Qr. Then since DaU, a — 2, • • • , d, satisfies 
V{DaU) = in Qr, by the above lemma, it follows that 

II^-C)oHIl2(Q.i) < ^ll-C'aHlL2(Q.2) < ^IIH1l2(Q«), 

where ri < r2 < R and N depends only on d, m, 5, and radii ri, r2, R. If we 
further consider Dapu, a, (3 = 2, ■ ■ ■ ,d, which satisfies V{Dapu) = in Qr, again 
by the above lemma 

\\DDapu\\L^(^Q^^) < N{d,m,S,ro,ri)\\Dai3u\\L^(^Q^^), 



8 



H. DONG AND D. KIM 



where ro < ri- By combining the above two incquaUties we obtain 
\\DD^0u\\L,iQ^^) < N\\Du\\L^^Q^y 

where we may also have || itlUaCQfi) instead of ||£'it||i2(Q^) on the right-hand side. 
By repeating the same reasoning, we have 

\\DD',,u\\L,iQ^) < N\\Du\\L,iQ,), 

where A: is a positive integer, r < R and N = N{d,m,S,r,R,k). Considering the 
derivatives of u in time as weh, in general we have the following lemma. With an 
additional assumption that the coefficient matrices are symmetric, a similar result 
was proved in [p^ . 

Lemma 3.3. Let < r < R < oo and A"^ = A"'^{xi), a, P ^ 1,2, - ■ ■ ,d. Assume 
u £ Cf^^ satisfies 

ru=0. (3.2) 

in Qr. Then we have 

\\D\D'i,u\\L,(Q^) + \\D\Di,DML.i^Q.) < N\\Du\\l,^q^), 

wherei,j are nonnegative integers satisfying i + j > 1 and N — N{d,m,6,R,r,i,j). 

Proof. Note that 'P{D\D-'^,u) — in Qn. Thus, thanks to the argument shown 
before this lemma, it suffices to prove 

||%||l.(q.) < N{d,m,5,R,r)\\Du\\L,(Q^). 

Set 

" R-r 

TQ^r, rn = r + ^ , n = l,2,---, 

fe=i 

We choose (^n{t,x) G such that 

1 on Q„ 



and 



on M'*+i\(-4,4) xi3,„ 
|-DC»I < ^- 



R-r 
Also set 

An = KIL,(Q(")), B = \\Du\\l,(Q^). 

After multiplying both sides of (3.2) from the left by (^{u^, ■ ■ ■ ,m™) and inte- 
grating on Qji, we get 

d ^ dm-. 

{u\^nfdxdt^ D^{ulC)AffD0u'dxdt = O. 

i=l •'Qr a,l3=l ij = l •'Qr 

From this and the Young's inequality, we have 

V / (v},Ufdxdt 

d in ^ 

Y (Da.ulA^fDt,u^C + 2ulA^fD0U^CnD^(:n)dxdt 



a., (3—1 — 1 
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'[ \DtLt\\Du\dxdt+l-y2 [ {ulC„fdxdt + N2^''[ \Du\^dxdt, 

where N — N(d, m, 5, R, r). This inequahty yields, for any e > 0, 
An <e||i^iXt||^^(Q(„,)+7V(d,TO,5,i?,r)(2" + e-i)B. 



(3.3) 

Since % also satisfies (3.2) in Q^j, by Lemma p.2[ , we have 

||i?t^lL,(Q(„,) < Afi2"||%|L,(Q(„+i,) = iVi2"A„+i, (3.4) 

where iVi — Ni{d,m, 6, R,r). Upon setting e = l/(3A^i2"), we get from (3.3) and 
(0) that 

A„ < A„+i/3 + iV2"B. (3.5) 
We multiply both sides of (|3.5|) by 3^" and sum over n to obtain 



n=0 



n=0 



Therefore, 



^3-"A„ < ^3-""iA„+i +7V^(2/3)"B. 

oo 

Ao < 7V^(2/3)"B. 



n=0 



The lemma is proved. 



□ 



Owing to the structure of the divergence form systems, the same type of inequal- 
ity as in the above lemma holds true if u in the left-hand side is replaced by U, the 
definition of which is 



d m 



[/ ^ Ai^T^^tt, i.e., = i = l,---,m. 

(3=1 /3=lj = l 



Lemma 3.4. Let = A°''^{xi). Assume u e C^^ satisfies (|^ in Q^. Then, 
for nonnegative integers i,j 

PP^C^IIl.(q.) + \\DlDi,D,U\\L,iQ,) < ^II^HIl.(Q4), (3-6) 
where N = N{d, m, d, i, j) > 0. 

Proof. As before, to prove ( |3.6D it is enough to show 

\\U\\l,(q,) + \\DiU\\l,^q,) < N\\Du\\L,^Q,y (3.7) 



Indeed, if this holds true, due to the fact that DlD-!^,u also satisfies (3.2) we have 

\\DlDi,UU,^Q,^ + \\DlDi,D,U\\L,iQ.) < N\\DlDi,DuU,^Q,y 

Then by Lemma 3^ we bound the right-hand side of the above inequality by a 
constant tim es | |-DmHltCO^'Ii so we arrive at the inequality (3.6). 
To prove (3/7), we observe that in Q4, 

d d 
Q=2,/3=l a=2,l3=l 

where the last equality is due to the independency of A"-^ in x' G M"^^^. Therefore, 

\\DiU\\l,(q,) < N\\ut\\L,iQ,)+N\\D,^,u\\L^^Q^y 
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By Lemma [3?^ , the right-hand side of the above inequahty is bounded by a constant 
times \\Du\\L^(^Q^y Thus 

\\DMl.AQ.) < N\\Du\\L,^Q^y 

It is clear that 

\\U\\L,iQ,) < N\\Du\\L,^Q,y 
Therefore, the inequahty ( |3.7j) and thus Lemma ^.4| are proved. □ 

As usual, for /i £ (0, 1) and a function w defined on I? C K'^+^, we denote 

\w{t,x) - w{s,y)\ 



Lemma 3.5. Let A"'^ = A"^(a:i). Assume that u e Cf^^ satisfies ( P^ m (54. 
T/ien we have 

[U]cr;2^Q,) < mDuh.iQ.), (3.8) 
[i5.'M]ci/=(Q,) <^II^"IU2(Q4): (3-9) 
where N = N{d, m,5) > 0. 

Proof. We first prove ( |3.8| ). By the triangle inequality, we have 



(i,a;),(s,a)eQi I* - S\^'^ + \x - y\^/'^ ~ a:^,y^e{-l^),'x^^y^ \xi - Vl]^/'^ 

\U(t,y,,x')-U{s,yi,y')\ _ ^ , ^ 

(t,a:'),(s,!/')6Qi 
(t,a;')5^(s,;/') 

Hence the inequality (3^) follows if we prove li < \\Du\\]^^(^Q^y i — 1,2. 

Estimate of Ii: By the Sobolev embedding theorem U(t,xi,x'), as a function of 
xi G (—1, 1), satisfies 

\U{t,xi,x')-U{t,yi,x')\ ^,1,^,^^ .„^„. 
sup ^ -Yj^ <N\\U{t,-,x)\\wU-i.i)- (3-10) 

On the other hand, there exists a positive integer k such that U{t, xi, x') and 
DiU{t,xi,x'), as functions of {t,x') £ Q[, satisfy 

sup {\U{t,xi,x')\ + \DiU{t,x,,x')\) 

< l|C^(->a;i,-)||vi/^(Qi) + \\DiU{-,xi,-)\\wi(Q[)- 
This implies that, for all {t,x') £ Q[, 

\U{t,xi,x')\'^ dxi + [ \DiU{t,xi,x')\'^ dxi 

lL(Q2)- 



i<ijl+j2<k 

This combined with ( [3.10| ) shows that 

h<N \\DlDl^Di^,UU,iQ,j < N\\DuU,^Q,), 

*<1 J'l+J2<fc 
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where the last inequaUty is due to Lemma 3.4 



Estimate of I2: Again using the Sobolev embedding theorem, we find a positive 
integer k such that U{t,yi,x'), as a function of {t,x') e Q[, satisfies 

\U{t,yi,x')-U{s,yi,y')\ 

+ ^mUi;y.,-)\\^,,,,y (3.11) 

For each i,j such that i+ j < k, D\D-'^,U{t,yi,x'), as a function of j/i e (—1,1), 
satisfies 

sup \D\Di,U{t,y^,x')\ 
yie(-i,i) 

<A^||I?jl?^,,C/(t,.,xOIU,(-ia)+A^I|i^ji^'.'^if/(i,-,a:')IU.(-i,i)- 
This together with ( 3.1l| ) gives 



where the last inequality follows from Lemma 3.4. Hence the inequality (3. 



proved. The proof of (|3.9D is done by repeating the same reasoning as above with 
the help of Lemma |3.3| . □ 

By using a scaling argument, wc have the following corollary. 

Corollary 3.6. Let r e (0, 00), k e [4, 00), and A"''^ = A°'f^{xi). Assume u e C^^ 
satisfies (|3.2|) in Q^r- Then we have 



{\U~{U)QA)Q^<Nn-'/^{\Du\')]/l 



(3.12) 
(3.13) 



where N = N{d,m,d) > 0. 



Proof. We prove only ( 3.12 ) . The inequality ( 3.13 ) is proved similarly. By a scaling 
argument, i.e., by considering u{r^t,rx) and A°'^{rxi), it suffices to prove (3.12) 
when r = 1. In this case, to use again the same type of scaling argument we define 

u{t, x) = u{{k/4)\ iK/A)x), A^l^ixi) = A"'5((K/4)a;i). 

Since Vu ~ in Q„, we have 

-ut + DaiA^^Dpii) = 

in Q4. Then by Lemma |3.5| applied to IJ — X]/3=i A^'^Di^u, we have 

[U]c^/HQ,)<N\\Du\\L,iQ,)- 

Note that 

{\U-iU)QA)Q, < [U]cu^Q,) < [f^]cv.(Q,,,) 

The corollary is proved. □ 
The following is the main result of this section. 
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Proposition 3.7. Let k G [8,oo), r e (0,oo), A"'^ = A°'P{xi), and g E L2,ioc- 
Assume that u G joc satisfies 

Vu = div g 

in Q^r- Then we have 

where N depends only on d, m, and 5. 

Proof. By performing the standard mollifications, we may assume it, g and 
are smooth. Take C G Co° such that 

Q — 1 on QKr/2, C = outside (— (^r)^, (^r)^) x B^r- 



By Lemma 3.1 (iii), there exists a unique w e ?i2((— (kt")^, 0) x M'') satisfies 

Vw ^ div(Cg). 

Moreover, 

\\Dw\\L2{Q^r) < IIC5llL2((-(Kr)2,0)xR'') < I|5|Il2(Q..)- (3-15) 

Let V = u — w so that Vv = in QK.r/2- Note that by the classical result w is in 
fact infinitely differentiable in Q^r because the coefficients of the operator as well 
as Cg are smooth. Hence v is also infin itely differentiable in Q^r- 
Denote V = A^I^Dpv. By Corollary |^ we have 

(1^ - (^)qJ)q„ + {\D,'V-{D^,v)Q^)Q^ < Nn~'/^{\Dv\')lll. (3.16) 
Now we observe that 

i\U-iU)Q.\)Q^ + m'U^ {D.,,u)qJ)q^ 

<2{\U- {V)q,.\)q^ + 2 {\D,>u~ {D.,'v)q^\)q^^ . 

Therefore, upon using ( 3.15 ), ( 3.16 ) and the triangle's inequality, we bound the 
left-hand side of ( |3.14D by 

2{\V-{V)Q,.\)Q^^+2i\D,,v- {D,,v)q,.\)q^+N{\Dw\)q,. 

<iV.-V2(|^^|2)V2^^^^(.+2)/2(|^|2)V2^ 
<iV.-V2(|^,|2)l/2^+^^(.+2)/2(|^|2)V2^^ 

The proposition is proved. □ 

We will also make use of a generalization of the Fefferman-Stein Theorem proved 
recently in |Q. Let C„ = {C'„(io, «i, • • • , id), io, - ■ ■ id ^ Z}, n G Z be the filtration 
of partitions given by parabolic dyadic cubes, where 

Cniio,ii, ■ ■■ ,id) 

= [zo2-2", (lo + 1)2-2") ^ [^^2-", (ii + 1)2-") X ... X [i,2-", (z^ + 1)2-"). 
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Theorem 3.8 (Theorem 2.7 [||). Let p e (0, 1), F,G,H e Li. Assume G > \F\, 
H > and for any n E Z, and C E Cn there exists a measurable function 
defined on C such that |i^| < F^ < G on C and 



j \F^ - {F'^)c\dxdt < [ Hdxdt. 
Jc Jc 



Then we have 

\\FrL^<N\\HU^\\G\\l-\ 

provided that H, G € Lp. 



4. Proof of Theorem 2.2 



In this section we complete the proof of Theorem 2.2 



Lemma 4.1. Let k>8, r>0, AeA, tl^E"^,uE Cf^^ and g G L2joc- Assume 

- nt{t,x) + Dk (^A''^{y)Diu{t,x)'^ ^ div g, (4.1) 
where y — ipix), 4> — 4'^'^, and 

d 

A'^\y)= Dy^<j,u{y)A'-P{y,)Dy^cj,i{y), k,l = \,---,d. 
Then there exist constants v = v{d, 5) >1 and N = N{d, m, S) > such that 

d 

{\ju - (J[/)qJ)q^ + - iJ^M)Q^ 

0=2 

where 

up{t,x)^{Dy^v){t,xP{x)), v{t,y)^u{t,4>{y)), (4.2) 

d 

J{y)=det{dcp/dy), U{t,x) = YA^''{M^))Mt,^)- (4-3) 

0=1 

Proof. Without loss of generality, we assume ■0(0) = 0. From the integral formula- 
tion of ( [4.1| ) , we see that v satisfies 

~{Jv)t+Dy^ (A''^{yi)JDy^,v) = div5, 

where 

d 

13=1 

So Jv satisfies 

-{Jv)t + Dy^ (A'-('{yi)Dy^{Jv)) -divg, 

where 

9= (9i,-- - 

d d 

ffo = ^E5/3(i>'^(y))P^^«)Wy)) + E^"^(yl)^(^'y)^^'.•^■ 
/3=l 13=1 
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Then by Proposition 3.7, 

< Nn-'/'{\D{jv)n]^i+Nn^''+'yHmiil, (4.4) 

where 



y(<,y) = ^ii^(yi)i?,,(Ji;). 

/3=1 

Now we observe that 

d 

{\ju - {ju)qJ)q^ + i\Ju0 - {Jy'P)Q.\)c 
0=2 



< 2 {\JU - {V)q^A)q^ +2^ {\Jup - (D,,(J^))q„J)q^ := 2(/i + h), 



f)=2 



where > 1 is a constant obtained in the foUowing observation. There exist 
constants v as weh as N depending only on d and 5 such that, for a nonncgative 
measurable function f(t,x), 



f{t,x)dxdt<N -f f{t,(j){y))dydt, 
Qp JQ.p 

fit, 0(y)) dydt<N -f fit, x) dx dt. 



Thus 



h = 



\Ji^Pix))Y,A"'iM^))iDy,v)it,^l,ix)) - iV)Q^Jdxdt 



p d 

<N f \J{y)Y,A'^iyi)iDy,v)it,y) - iV)Q,Jdydt 
< N i\V - iV)Q^^)^^+ N{\v\)q^^. 



Similarly, 



h= f \Ji'^{x))iDy^v)it,^ix)) ~ iDy^iJv))Qjdxdt 



<N { \Jiy)iDy^v)it,y) - {Dy^{Jv))Q^^\dydt 



< 



Ni\Dy.iJv) ~ {Dy>iJv))Q^A)^^^+Ni\v\)Q^^. 



(4.5) 



Using the above two sets of inequalities for Ii and I2 as well as using (4.4) with vr 
in place of r and (4.5), we obtain the desired inequality in the lemma with v'^ in 
place of V. The proof is completed upon simply replacing by another constant 

V. □ 



With the aid of Lemma 4.1, we estimate the mean oscillations of JU and J up 
for general operators. 
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Proposition 4.2. Let 7 > and r, cr e (l,oo) satisfy 1/t + 1/ct = 1. Let v — 
i^{d,S) > 1 be the constant in Lemma \4-^ , B" = B" = C = 0, and g G L2.ioc- 
Assume that u G vanishes outside Qji for some R G (0, Rq] and satisfies 

Vu = diY g. Then under Assumption 2.1 ("/), for each r £ (0,oo), k > 8, and 
{to,xo) E R'^^"'^, there exist a diffeomorphism ij] £ , coefficients A^^,(3 = 1, - • • ,d 
(independent of u), and a positive constant N — N(d, m, S, r) such that 

d 

{\JU - (^C^)Q.(to,=.o)l)Q,(to,,„) + - (•^^/3)Q.(*o,xo)l)Q^(t„,,„) 

13=2 



where Up, J, and U are defined as in (4.2) and ( |4.3| ). 

Proof. We fix a K > 8 and r G (0, 00). Choose Q to be Q^Krito, x^) if vnr < R and 
Qu if i/Kr > R. Let {t*,x*) be the center of Q and y* ~ ip{x*). By Assumption 
2j| (7), we can find ijj e and A = A{r) G A satisfying (|2^. We set 

d 

A''^{yi)= 5] i?feVc.(a:*)A'='(2;i)AV'/3(a^*), a,P = l,---,d, 

k,l = l 

and 

d 

A'^'iy)^ DyMyW^{yi)DyMy)^ k,i = ir--,d, 

a,/3=l 

where y = tpi^)- The ellipticity constants of A and A may not be S, but they 
depend only on 6. Note that 



^Ut + Dk{A^'Diu)=divg+Dk {A^'Diu- A^'Dm) . 



Thus by Lemma 4.1 with a shift of the coordinates, 

d 



{\JU-{JU)q 

13=2 

where v = v{d, 5) >1, N = N{d, m, S) > 0, and 



/ Ig^l^dxdt:^ \A''' Diu- A''' Diu\^ dx dt 

<2 f \A''^DiU- A''^Diu\^ dxdt + 2 f \A''^DiU- A'^'Diu]^ dxdt 



I QuKr{to.Xo) J Q,^Kr{to,Xo) 

:=2(/i+/2) 

Note that 



d 



A''\y^) = J2 DyM{y*W^{yi)DyMy*)- 

Q,/3 = l 

Thus by (O) 
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llioo(Q 



JQ 



\Du\ dxdt 



< NW j \Du\-'dxdt. 
On the other hand, by the Holder's inequahty, we have 

h < Nll{^lli\ 

where 

\A'^\Mx)) - A''\t,x)\^" dxdt, 



(4.8) 



(4.9) 



22 



\Du\'^'' dxdt. 



Due to Assumption 2.1 (7), 



/21 < N-f\Q\ < N{iyKrf+^j 



This together with (4.7)-(4.9) yields ([4. 61). The proposition is proved 



□ 



The next corollary follows immediately from Proposition 4.2 by using the triangle 
inequality. 

Corollary 4.3. Let •j > 0, k > 8 and t, a E {1,00) satisfy 1/t + l/a ~ 1. Suppose 
that B°' — B°' = C = 0; g E L2.10C, o,nd u € vanis hes outside Qn for some 
R G (0,i?o] satisfying Vu = div g. Under assumption \2.\ (^), for each n G Z 
and G G C„, there exist a diffeomorphism G coefficients A^^ , f3 = 1, - • ■ ,d 
(independent of u), and a constant N = N{d, m, 6, r) such that 



{\JU - {JU)c\)c + E (l'^"/' - iJ^p)c\)c < N{H)c, 

/3=2 



where Ufj, J and U are defined as in (4.2) and ( |4.3| ), and 



d + 2 

H = K— 



(4.10) 



Proposition 4.4. Let p G {2, 00). Assume B°' = B°' = C = 0. Then there exist 
positive constants 7, N and R G (0, 1] depending only on d, ni, p, and 5 such that 
under Assumption 2.1 ('y), for any u G vanishing outside Qrr„ and g G Lp, 
we have 

\\DuU,<N\\gU,+N\\uU,, (4.11) 

provided that Vu— div g. 

Proof. Let 7 > 0, k > 8 and R G (0, 1] be constants to be specified later. Let 
T = (p + 2)/4 > 1 such that p > 2r. We take n G Z, C G C„ and let G * be the 
diffeomorphism from Corollary 4.3 corresponding to the chosen n and C. We also 
obtain corresponding Up, J and f/ as in (4^) and (^). 
It is easily seen that 



\Du\ < N2{d,S)J2\Jui3\+ N2id,S)\JU\ < N3id,S)\Du\. 

0=2 



PARABOLIC AND ELLIPTIC SYSTEMS 



17 



We set 



F=\Dul ^N2^\Jup\ + N2\JU\, G^N:i\Du\. 

/3=2 

By the triangle inequality and ( 4.1(]| ), 

{\F^ -{F^)c\)c<N{H)c. 



where H is defined in Corollary 4.3. Now by Theorem 3.S, we get 

\\Du\\l^ = \\F\\l^ < NWHUJCrC < N{e)\\H\\l + e\\G\\l 
Upon taking a small e > 0, it holds that 

\\Duh,<N\\H\\L,. 



(4.12) 



We use the definition of H and the Har dy-Littlewood maximal function theorem 
(recall p > 2r > 2) to deduce from (|4.12|) 



\Du\\l^<Nk^ {\\g\\ 



Hkp) + N{k^j^ + K^RRo + K 



\Du\\l,. 
(4.13) 

By choosing n sufficiently large, then 7 and R sufficiently small in (4.13) such that 

N{K—-f^ + K—RRa + K"5) < 1/2, 
we come to ( |4.11 ). The proposition is proved. □ 



Proof of Theorem Thanks for the duality argument, it suffices to prove the 
case p > 2. For T — oo, the theorem follows from Proposition 4.4 by using a 
partition of unity and an idea by S. Agmon; see, for instance, the proof of Theorem 
1.4 1^^. For general T S (— cx),cx)], we use the fact that u = w for t < T, where 
w e Hp solves 

Vw — Xw = xtKriVu — Xu). 
This finishes the proof of the theorem. □ 



5. A REMARK ABOUT ELLIPTIC SYSTEMS 

For elliptic systems, the condition on diffeomorphisms ip and 4> can be relaxed. 
Indeed, we only require and to be in C*^'^ and Dip has locally small mean 
oscillations. More precisely, we impose the following assumption on ip and A, 
which is weaker than the one in Section |^. 

Let vj/ be the set of C°'^ diffeomorphisms -0 : M'' — > M'' such that the mappings 
Tp and (p — ip^^ satisfy 

\Dip\<5-^, \D(P\<8-^. 

Assumption 5.1 (7). There exists a positive constant Rq E (0, 1] such that, for 
any ball B of radius less than Rq, one can find anAEA and aip — {tpi , • • • , G ^' 
such that 

V/ \A''\M^))"A''''Do.ADf3iJiJ{x)\dx<-f\B\. (5.1) 
k,i •'B 
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Lemma 5.2. Let k > 8, r > 0, A e A, e , u e C;^c(M'') and g G L2Joc(K'^). 
Assume 

(A''\y^)Dy^4>^{y)Dy^4>p{y)J-^Dpu{x)) - div g, 

where y ~ ipi^) and (j> ~ ^p^'^ . Then there exist constants v — v{d,5) > 1 and 
N — N(d,m,d) such that 

d 
13=2 



where 



J{x) = det{dij/dx)-\ U{x) = A^'^(tlji{x))up(t,x) 



(5.2) 
(5.3) 



Proof. The proof is similar to that of Lemma 4.1. From the integral formulation, 
it is easy to see that v satisfies 

DyAA''^{y')Dy^,v) = Dy^ [JDpi^^gp) . 

The lemma then follows from Proposition 3/7. □ 

Proposition 5.3. Let 7 > and T,cr ^ (l,oo) satisfy 1/r + \/cr = 1. Let v = 
:^(d, (5) > \ he the constant in Lemma ^.1, i?" = _B" ~ C = 0, and g E L2joc{^'^)- 
Assume that u G C^{W^) vanishes outside Br for some R £ (0,i?o] and satisfies 
Cu = div g. Then under Assumption \5.\ (7 ), for each r G (0, 00), k > 8, and Xq G 
M.'^ , there exist a diffeomorphism ip E , coefficients A^^ , /3 = 1, • • • ,d (independent 
of u), and a positive constant N — N(d, m, (5, r) such that 

d 



(i=2 



|2\l/2 



|2\l/2 



^ ^^(d+2)/2^1/(2.) (i^^i 



2r\l/(2r) 



(5.4) 



where Up, J and U are defined as in (5.2) and ( |5.3| ). 

Proof. We fix a K > 8, and r G (0,oo). Choose i? to be i?yKr(^o, 2^0) if '^'f < R 
and Sfl if vKr > R. By Ass ump tion 5.1 (7), we can find ?/; G ^' and A = ^(s) G A 
satisfying (5J^). By Lemma |5.2| with a shift of the coordinates, 

d 

(|C^- (C^)b.(.o)I)b^(,„) +E ("/3)s.(^o)I)b,(,„) 



/3=2 



j2n1/2 



|2\l/2 



'B,^,.(xo) 



<Nn^''+^)/^\g+g\% 
where A'' = N{d, to, (5) > and 

9c. = {A'^\yi)Dy,<p^{y)Dy,(l>p{y)J-^ - A^^{x)) Dpu{x). 



(5.5) 
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By the definition of A, 

\gJ^dx<N I \{A''\i,i)^ A^''D^i:kDp4,iJ)\^\Du\^ dx 



where 



< iV/y^/y", (5.6) 



\Du\^^ dx 



Due to Assumption 5.1 (7), 



/i < V / lA'^Hi^i) - A"l'Da,ibkD0iJ;iJ\^''dx < N^\B\ < N{vnrf-f. 

, JB 



This together with (|5.5| ) and ( |5.6| ) yields (|5.4| ). The proposition is proved. □ 



FoUowing the arguments in the previous section, we obtain the result of Theorem 



2.5 under Assumption 5.1. We omit the details. 
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